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Abstract 

It is shown that exact, amphtude-based resummation allows IR-improvement of the usual 
DGLAP-CS theory. This results in a new set of kernels, parton distributions and attendant 
reduced cross sections, so that the QCD perturbative result for the respective hadron- 
hadron or lepton-hadron cross section is unchanged order-by-order in at large squared- 
momentum transfers. We compare these new objects with their usual counter-parts and 
illustrate the effects of the IR-improvement in some phenomenological cases of interest 
with an eye toward precision applications in LHC physics scenarios. 
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1 Introduction 



With the impending start of the LHC, the era of precision QCD at the LHC, by which 
we mean 1% or better precision tags on the theoretical predictions, obtains. In this new 
era for perturbative QCD in high energy colhding beam physics, we have the extremely 
challenging task of proving that a given theoretical precision tag does in fact hold to that 
level. All aspects of the standard formula for hadron-hadron scattering in perturbative 
QCD have to be examined for possible sources of uncertainty in the physical and technical 
precision components of any claimed total theoretical precision tag. In this connection, we 
recall that, in Ref. [1], we have found that the resummation of large infrared (IR) effects 
in the kernels of the usual DGLAP-CS [2-5] theory results in their improved behavior 
in the respective IR limit. This improved behavior, for example, results in kernels that 
are integrable in the IR limit and which therefore are more amenable to realization by 
MC methods (our ultimate goal) to arbitrary precision. In this paper, we review the 
results in Ref. [1] and we show that the new IR-improved DGLAP-CS theory leads to a 
new set of parton distributions and reduced cross sections, so that the exact predictions 
from QCD for hadron-hadron and lepton-hadron processes are unaltered order-by-order 
in perturbation theory. The advantage is better control on the accuracy of a given fixed- 
order calculation throughout the entire phase space of the respective physical process, 
especially when the prediction is given by MC methods. 

The paper is organized as follows. In the next section, we recapitulate the new IR- 
improved DGLAP-CS theory, as it is not very familiar. In Section 3, we illustrate the 
use of the new IR-improved theory with some examples of phenomenological interest. In 
Section 4, we address the issue of what it does to the standard calculational apparatus for 
perturbative QCD predictions of hadron-hadron and lepton-hadron scattering processes at 
large squared momentum transfer. Section 5 contains some concluding remarks. Review 
of exact, amplitude-based resummation theory is presented in the Appendix. 

2 IR-Improved DGLAP-CS Theory 

Specifically, the motivation for the improvement which we develop can be seen already in 
the basic results in Refs. [5] for the kernels that determine the evolution of the structure 
functions by the attendant DGLAP-CS evolution of the corresponding parton densities by 
the standard methodology. Here, we already stress that the attendant evolution equations, 
under Mellin tranformation, are entirely implied by those of the Callan-Symanzik-type [2] 
analyzed in Refs. [3,4] in their classic analysis of the deep inelastic scattering processes. 
Thus, henceforward, we shall refer to these equations as the DGLAP-CS equations. Con- 
sider the evolution of the non-singlet (NS) parton density function q^^{x), where x can 
be identified as Bjorken's variable as usual. The basic starting point of our analysis is the 
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infrared divergence in the kernel that determines this evolution: 



where the well-known result for the kernel Pqq{z) is, for z < 1, 

P,,{z) = (2) 

when we set t = Iu/^^/uq for some reference scale /iq with which we study evolution to 
the scale of interest /x. u Here, Cp = {N^ — l)/{2Nc) is the quark color representation's 
quadratic Casimir invariant where is the number of colors and so that it is just 3. This 
kernel has an unintegrable IR singularity at z = 1, which is the point of zero energy gluon 
emission and this is as it should be. The standard treatment of this very physical effect 
is to regularize it by the replacement 

(3) 



(1-z) (l-z) 

with the distribution ^-^J-^.^^ defined so that for any suitable test function f(z) we have 

A possible representation of 1/(1 — 2)+ is seen to be 
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[l-z 



^(l-e-^)+lne5(l-^) (5) 



with the understanding that e J, 0. We use the notation 9{x) for the step function from 
for a; < to 1 for a; > and 5{x) is Dirac's delta function. The final result for Pqq{z) is 
then obtained by imposing the physical requirement [5] that 

dzPqq{z) = 0, (6) 

which is satisfied by adding the effects of virtual corrections at z = 1 so that finally 

n.W = Ci^(7^^ + ^5(l-^))- (7) 



The smooth behavior in the original real emission result from the Feynman rules, 
with a divergent 1/(1 — 2;) behavior as 2; — > 1, has been replaced with a mathematical 

^Wc will generally follow Rcf. [6] and set /iQ ~ J^qcd without loss of content since dt ~ dt' when 
t = In/i^/Agj^^, t' = ln/i^//iQ for fixed values of Aqcd, ^J■o■ 
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artifact: the regime 1 — e < z < 1 now has no probabihty at all and at z = 1 we have a 
large negative integrable contribution so that we end-up finally with a finite (zero) value 
for the total integral of Pqq{z). This mathematical artifact is what we wish to improve 
here; for, in the precision studies of Z physics [7-9] at LEPl, it has been found that such 
mathematical artifacts can indeed impair the precision tag which one can achieve with 
a given fixed order of perturbation theory. An analogous case is now well-known in the 
theory of QCD higher order corrections, where the FNAL data on pt spectra clearly show 
the need for improvement of fixed-order results by resumming large logs associated with 
soft gluons [10-12] H For reference, note that at the LHC, 2 TeV partons are realistic so 
that z = 0.001 means ~ 2 — 3 GeV soft gluons, which are clearly above the LHC detector 
thresholds, in complete analogy with the situation at LEP where z = 0.001 meant ~ 100 
MeV photons which were also above the LEP detector thresholds - just as resummation 
was necessary to describe this view of the LEP data, so too we may argue it will be 
necessary to describe the LHC data on the corresponding view. And, more importantly, 
why should we have to set Pqq{z) to for 1 — e < z < 1 when it actually has its largest 
values in this very regime? 

By mathematical artifact we do not mean that there is an error in the computations 
that lead to it; indeed, it is well-known that this -I— function behavior is exactly what one 
gets at 0{as) for the bremsstrahlung process. The artifact is that the behavior of the 
differential spectrum of the process for 2 — s> 1 in 0{as) is unintegrable and has to be 
cut-off and thus this spectrum is only poorly represented by the 0{as) calculation; for, 
the resummation of the large soft higher order effects as we present below changes the 
z — >■ 1 behavior non-trivially, as from our resummation we will find that the ^z^-behavior 
is modified to (1 — 2)'^"^, 7 > 0. This is a testable effect, as we have seen in its QED 
analogs in Z physics at LEPl [7-9]: if the experimentalist measures the cross section 
for bremsstrahlung for gluons (photons) down to energy fraction eo, eo > 0, in our new 
resummed theory presented below, the result will approach a finite value from below as 
— eQ whereas the 0{as) -I— function prediction would increase without limit as — In eo. The 
exponentiated result has been verified by the data at LEPl. 

The important point is that the traditional resummations in N-moment space for the 
DGLAP-CS kernels address only the short-distance contributions to their higher order 
corrections. The deep question we deal with in this paper concerns, then, how much of the 
complete soft limit of the DGLAP-CS kernels is contained in the anomalous dimensions of 
the leading twist operators in Wilson's expansion, an expansion which resides on the very 
tip of the light-cone? Are all of the effects of the very soft gluon emission, involving, as 
they most certainly do, arbitrarily long wavelength quanta, representable by the physics 
at the tip of the light-cone? The Heisenberg uncertainty principle surely tells us that 
answer can not be affirmative. In this paper, we calculate these long-wavelength gluon 

^Notc that we do not address px resummation effects herein; we use this example to illustrate the 
point that, whenever the phase space for radiation is restricted in a given variable, soft gluon radiation 
will result in large logs that need to be resummed, altough the methods to make this resummation may 
vary in detail. 
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(a) (b) 
Figure 1: In (a), we show the usual process q q{l — z) + G{z); in (b), we show its 
multiple gluon improvement q q{l — z) + 6*1(^1) + ■ ■ ■ + z = ^j. 



effects on the DGLAP-CS kernels that are not included (see the discussion below) in the 
standard treatment of Wilson's expansion. We therefore do not contradict the results of 
the large N-moment space resummations such as that presented in Ref. [13] nor do we 
contradict the renormalon chain- type resummation as done in Ref. [14]. 

We employ the exact re- arrangement of the Feynman series for QCD as it has been 
shown in Ref. [15,16]. For completeness, as this QCD exponentiation theory is not gen- 
erally familiar, we reproduce its essential aspects in our Appendix. The idea is to sum 
up the leading IR terms in the corrections to Pgq with the goal that they will render 
integrable the IR singularity that we have in its lowest order form. This will remove the 
need for mathematical artifacts and exhibit more accurately the true predictions of the 
full QCD theory in terms of fully physical results. 

We apply the QCD exponentiation master formula in eq. (l80|) in our Appendix (see 
also Ref. [15]), following the analogous discussion then for QED in Refs. [8,9], to the 
gluon emission transition that corresponds to Pqq{z), i.e., to the squared amplitude for 
q — > q{z) + G{1 — z) so that in the Appendix one replaces everywhere the squared ampli- 
tudes for the Q'Q Q"'Q" processes with those for the former one plus its nG analoga 
with the attendant changes in the phase space and kinematics dictated by the standard 
methods; this implies that in eq.(53) of the first paper in Ref. [5] we have from eq. (1801) 
the replacement ( see Fig. [2] ) 
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where A = q, B = G, C = q and Va->b+c is the lowest order amphtude for q 
G{z) + g(l — z), so that we get the un-normahzed exponentiated result 



where [8,9,15,16] 



P,,{z) = CFFyFsiiy^'^^^i^ - (9) 



as 4:Cf , . 

Iq = Cp—t = —— (10) 







6. = ^ + ^(^-h fill 



and 



Here, 



" 2 ' TT ' 3 2 



2 

/?o = 11 - 2?^/ 
, where nj is the number of active quark flavors, 

Ce = .5772 . . . 

is Euler's constant and T{w) is Euler's gamma function. The function Fyfs{z) was already 
introduced by Yennie, Frautschi and Suura [17] in their analysis of the IR behavior of QED. 
We see immediately that the exponentiation has removed the unintegrable IR divergence 
at 2 = 1. For reference, we note that we have in Q) resummed the term^ Ofln'^fl — 
z)t^ag), n > £ > k, which originate in the IR regime and which exponentiate. The 
important point is that we have not dropped outright the terms that do not exponentiate 
but have organized them into the residuals /3m in the analog of eq.( l80|) . The application 
of eq. (l80!) to obtain eq.([9]) proceeds as follows. First, the exponent in the exponential 
factor in front of the expression on the RHS of eq. (l80|) is readily seen to be from eq. (!77|) . 
using the well-known results for the respective real and virtual infrared functions from 
Refs. [15,16], 

SUMjRiQCD) = 2asReBQCD + 2a,5QCD(/W) 

^^2C7^^tln^ + C^^t ■ ^^^^ 



2 V TT E 2n vr ' 3 2- 



^Following the standard LEP Yellow Book [7] convention, we do not include the order of the first 
nonzero term in counting the order of its higher order corrections. 
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where on the RHS of the last result we have already applied the DGLAP-CS synthe- 
sization procedure in the third paper in Ref. [15] to remove the coUinear singularities, 
InAg^^/m^ — 1, in accordance with the standard QCD factorization theorems [18]. This 
means that, identifying the LHS of eq. flHOj) as the sum over final states and average over 
initial states of the respective process divided by the incident flux and replacing that 
incident flux by the respective initial state density according to the standard methods for 
the process q — > g(l — z) + G{z), occurring in the context of a hard scattering at scale Q 
as it is for eq.(53) in the first paper in Ref. [5], the soft gluon effects for energy fraction 

< z = Kmax/E give the result, from eq. flHOl) . that, working through to the /^-level and 
using q2 to represent the momentum conservation via the other degrees of freedom for the 
attendant hard process, 

2tt J J (27r)4 

^^^1 ^ (l. \ f ^'^y ^fa-(pi-P2-fci)+/'°<^'"°" 4^5QCD(fc)[e-'^-''-l]} 

h ' J {2tiY 

' ' J 

^ pSUMiH,(QCD)(z) fr7> H J^Jiy{E^-E2)+f<''"^-- ^SQcu{k)[e-'y''-l]} 

J h (27r) 



1 d^P2 



(14) 



where we set Ei = i = 1,2 and the real infrared function S'qcd(^) is well-known as 
well: 

C fl\ - ^^^P f Pi , 

•JQCDI/CJ — I T T — I I DGLAP-CS synthesized UOJ 



and we indicate as above that the DGLAP-CS synthesization procedure in Refs. [15] 
is to be applied to its evaluation to remove its collinear singularities; we are using the 
kinematics of the first paper in Ref. [5] in their computation of Pba{z) in their eq.(53), 
so that the relevant value of is indeed Q^. It means that the computation can also be 
seen to correspond to computing the IR function for the standard t-channel kinematics 
and taking | of the result to match the single line emission in Pcq. The two important 
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integrals needed in ( fT^ were already studied in Ref. [17]: 



oo 



zE 1. 



(16) 



Iyfs{zE,0) 



'zE - ki 

When we introduce the results in f|T6l) into f|H|) we can identify the factor 

/ + / **'*'^'(*''(;^)"^'^) S - / + 

(17) 

where is the unexponentiated result in the first line of ([8]). This leads us finally to 
the exponentiated result in the second line of (IHl) by elementary differentiation: 



Pba = P'BAZ^''FYFs{l,)e-^ 



2^1 



Let us stress the following. In this paper, we have retained for pedagogical reasons 
the dominant terms in the resummation which we use for the kernels. The result in the 
first line of (IT^ is exact and can be used to include all higher order resummation effects 
systematically as desired. Moreover, we have taken a one-loop representation of for 
illustration and have set it to a fixed- value on the RHS of (fT^ . so that, thereby, we are 
dropping further possible subleading higher order effects, again for reasons of pedagogy. 
It is straight forward to include these effects as well [19]. 

Here, we also may note how one can see that the terms we exponentiate are not 
included in the standard treatment of Wilson's expansion: From the standard method [20], 
the N-th moment of the invariants T^^^ , z = L, 2, 3, i = q, G, of the forward Compton 
amplitude in DIS is projected by 

^q{^^l . . . qi^N] qn 

T^N = 

[ N\ ■ ■ ■ ^p^'^ 

where Xsj = <5^/(2gp) in the standard DIS notation; this projects the coefficient of 
1/ {2xBj)^ ■ For the dominant terms which we resum here, the characteristic behavior 
would correspond formally to 7q-dependent anomalous dimensions associated with the 
respective coefficient whereas by definition Wilson's expansion does not contain such. In 
more phenomenologically familiar language, it is well-known that the parton model used 
in this paper to calculate the large distance effects that improve the kernels contains 
such effects whereas Wilson's expansion does not: for example, the parton model can 



(19) 
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be used for Drell-Yan processes, Wilson's expansion can not. Similarly, any Wilson- 
expansion guided procedure used to infer the kernels via inverse Mellin transformation, 
by calculating the coefficient of (l/z)"' in Wilson's expansion, will necessarily omit the 
dominant IR terms which we resum. Here, we stress that we refer to the properties of the 
expansion of the invariant functions Tj, not to the expansion of the kernels themselves, as 
the latter are related to the respective anomalous dimension matrix elements by inverse 
Mellin transformations. We need to emphasize that we are not saying that there is an 
error in the standard use of Wilson's expansion. We are saying that the IR improvement 
which we exhibit is not a property of that expansion as it is conventionally defined. 

The normalization condition in eq.Q then gives us the final expression 



CfFyfs{1c, 



(1 - z)> - /,(7,)5(1 - z) 



(20) 



where 



fqilq) 



2 1 
+ 



The latter result is then our IR-improved kernel for NS DGLAP-CS evolution in QCD. 
We note that the appearance of the integrable function (1 — z)"^'^^'^ in the place of ^^^J^^^^ 
was already anticipated by Gribov and Lipatov in Refs. [5]. Here, we have calculated the 
value of 7g in a systematic rearrangement of the QCD perturbation theory that allows 
one to work to any exact order in the theory without dropping any part of the theory's 
perturbation series. 

The standard DGLAP-CS theory tells us that the kernel PGq{z) is related to Pqq{l — z) 
directly: for 2; < 1, we have 



Gq\ 



CpFvFsilqy 



2"1- 



1 + (1 



(22) 



This then brings us to our first non-trivial check of the new IR-improved theory; for, the 
conservation of momentum tells us that 



/ dzz{PGq{z)+P^^{z))=Q. 

Jo 



(23) 



Using the new results in eqs. ( l20|22l) . we have to check that the following integral vanishes: 

dzZ i _Li_^;27, + _ ^)7, _ /^(^^)5(1 -z)] (24) 



To see that it does, note that 



1-z 



1+1 ^ 1 
= -1 + ^ — 



1-z 



(25) 
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Introducing this result into eq. (l2^ we get 



Jo i- — z 

The integrals over the first two terms on the right-hand side (RHS) of fl26l) exactly cancel 
as one sees by using the change of variable 2; —> 1 — 2; in one of them and the integral over 
the last two terms on the RHS of vanishes from the normalization in eq.([2]). Thus 
we conclude that 

/ = 0. (27) 

The quark momentum sum rule is indeed satisfied. 

Having improved the IR divergence properties of Pqq{z) and PQq{z), we now turn to 
Pgg{z) and PqG^z). We first note that the standard formula for PqG{z), 

PqG{z) = \{z' + {l-zfl (28) 

is already well-behaved (integrable) in the IR regime. Thus, we do not need to improve it 
here to make it integrable and we note that the singular contributions in the other kernels 
are expected to dominate the evolution effects in any case. We do not exclude improving 
it for the best precision [19] and we return to this point presently. 

This brings us then to Pgg{z)- Its lowest order form is 

Pgg{z) = 2Cg{^ + + z(1 - z)) (29) 

z 1 — z 

which again exhibits unintegrable IR singularities at both z = 1 and z = 0. (Here, Cg 
is the gluon quadratic Casimir invariant, so that it is just = 3.) If we repeat the 
QCD exponentiation calculation carried-out above by using the color representation for 
the gluon rather than that for the quarks, i.e., if we apply the exponentiation analysis in 
the Appendix to the squared amplitude for the process G G{z) + G(l — z), we get the 
exponentiated un-normalized result 

Pgg{z) = 2CGFYFs{lG)ei'^ f + -^(1 - z)^« + Uz^+^'^il - z) + z{l - zY+^o)\ 

\ z 1 — z 2 / 

(30) 

wherein we obtain the 7g and 5g from the expressions for 7^ and 5q by the substitution 
Cf Cg'- 

IG = CG^t = ^ (31) 

We see again that exponentiation has again made the singularities at z = 1 and z = 
integrable. 



Cg 


—t -- 

n 


^Cg 


IG 


a 

+ - 


sCg.71^ 


2 


71 ^ 3 
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To normalize Pgg? we take into account the virtual corrections such that the gluon 
momentum sum rule 



[ dzz{2nfPgG{z)+PGG{z))=0 
Jo 



(33) 



is satisfied. This gives us finally the IR-improved result 

Pgg{z) = 2CGFYFs{lG)e^'n^^z^'' + - ^)^^ 

+ \{z^+^^{l -z)+ z{l - zY+^^) - /g(7g)5(1 - z)} 

where for fcilG) we get 

1 f=-2^G 



(34) 



/g(7g) = r. r . ^""^ + V77^ ^ + 



QCgFyfs{ig) 7g(1 + 7g)(2 + 7g) (1 + 7g)(2 + 7g) 

1 1 

^ 2(3 + 7g)(4 + 7g) ^ (2 + 7g)(3 + 7g)(4 + 7g) ' 



(35) 



It is these improved results in eqs.( !20]l22p4l) for Pqq{z), PGq{z) and Pgg{z) that we use 
together with the standard result in (l29l) for PqG{z) as the IR-improved DGLAP-CS 
theory. 



For clarity we summarize at this point the new IR-improved kernel set as follows 

'l + z^ 



1- z 



.(l-^)>-/,(7,)5(l-^) 



(36) 



P^7(z) = C^Fy^,(^,)ei^^i±ii^z^% (37) 

P^^J(^) = 2CGFyMlG)e'''H^^^''' + - 

+ liz'^'H^ - ^) + ^(1 - zY^'n - /g(7g)5(1 - z)}, (38) 

P,g{z) = liz' + {1 - zn (39) 

where we have introduced the superscript exp to denote the exponentiated results hence- 
forth. 

Returning now to the improvement of PqG{z), let us record it as well for the sake of 
completeness and of providing better precision. Applying eq. (l80|) to the process G Q+Q, 
we get the exponentiated result 

P;^^iz) = FyM^G)e'^'o^{z'il - Zyo + (1 _ ^)2^7C|. (40) 

The gluon momentum sum rule then gives the new normalization constant for the P^q 
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via the result 



Uf 1 



•^''^^''^ Cg(1 + 7g)(2 + 7g)(3 + 7g) ^7g(1 + 7g)(2 + 7g) ^ (1+7g)(2 + 7g) ^^^^ 



1 1 



2(3 + 7g) (4 + 7g) (2 + 7g) (3 + 7g) (4 + 7g) ' 

The constant fc should be substituted for fc in P^q whenever the exponentiated result 
in (HOl) is used. These results fl39l) . fHOjl . and pij) are our new improved DGLAP-CS 
kernel set, with the option exponentiating Pgc as well. 

In the discussion so far, we have used the lowest order DGLAP-CS kernel set to 
illustrate how important the resummation which we present here can be. In the litera- 
ture [21,22], there are now exact results up to 0{al) for the DGLAP-CS kernels. The 
question naturally arises as to the relationship of our work to these fixed-order exact re- 
sults. We stress first that we are presenting an improvement of the fixed-order results such 
that the singular pieces of the any exact fixed-order result, i.e., the ^-^J-^^^ parts, are ex- 
ponentiated so that they are replaced with integrable functions proportional to (1 — z)'^~^ 
with 7 positive as we have illustrated above. Since the series of logs which we resum 
to accomplish this has the structure Q;ft^ln"(l — z), i > n these terms are not already 
present in the results in Refs. [21,22]. As we use the formula in eq. (!80|) . there will be 
no double counting if we implement our IR- improvement of the exact fixed-order results 
in Refs. [21,22]. The detailed discussion of the application of our theory to the results 
in Refs. [21,22] will appear elsewhere [19]. For reference, we note that the higher order 
kernel corrections in Refs. [21,22] are perturbatively related to the leading order kernels, 
so one can expect that the size of the exponentiation effects illustrated above and below 
will only be perturbatively modified by the higher order kernel corrections, leaving the 
same qualitative behavior in general. 

In the interest of specificness, let us illustrate the IR-improvement of Pgg when cal- 
culated to three loops using the results in Refs. [21,22]. Considering the non-singlet case 
for definiteness (a similar analysis holds for the singlet case) we write in the notation of 
the latter references 

oo 

p+ ^ pv , pv — \^ ( ^\n+l p{n)+ /^2^ 
ns qq ' qq / j \ a J ns \ J 



"-An' 

=0 



where at order 0{as) we have 



pW+(^) = 2C^{-i±^ + 1^(1 - z)} (43) 
{l-z)+ 2 

which shows that PnJ~^{z) agrees with the unexponentiated result in (I7j) for Pgg except 
for an overall factor of 2. We use this latter identification to connect our work with that 
in Refs. [21,22] in the standard methodology. In Refs. [21,22], exact results are given for 
Pnl^~^{z), and in Refs. [22] exact results are given for Pj^^^(z). When we apply the result 
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in (IHOl) to the squared amplitudes for the processes q^q + X,q^q + X', we get the 
exponentiated result 

a (44) 

where Pgg^{z) is given in (!39l) and the resummed residuals Pns~^ , i = 1,2 are related to 
the exact results for Pn]^ , i = 1, 2, as follows: 

(^) = ^if(^) - Br^Al -z)+ A^l^z) (45) 



where 



A«+(.) = -4C^7r5ai±^ - - z)} 

Ai^i+W = -4CH-5i)^{^ - M(l - .)} ^^^^ 



and 



(47) 



52 = fi2 + 4C^vr5i/, 

53 = ^3 + 4C^(7r5i)V, - 27r5i52. 

Here, the constants Bi, i = 2,3 are given by the results in Refs. [21,22] as 

B, = 4CgCf{^ + yC2 - 3C3) - 4C^n;(l + ^(2) + 4C|(^ - 8(2 + 6(3) 

^ ^ -^M 54^ 20^^ 18^^ 
+ IQCaCU'-^ + C2C3 - - ^C2^ + ^Cs + yC5) 

+ 16C^CH-^ + ^C2-^C2^-^C3 + ^C5) (48) 
^ 576 27 8^ 9 2^ 

144 27 9 ^ 

^ / 23 5, 29,0 X7^> 

+ IQClnpi + — C2 + — Cl Ca 

^ ^ 16 12^ 30^^ 6 ^ ' 

+ 16C|(| - 2C2C3 + ^C2 + fa + ^Cs - 15C5), 

where C„ is the Riemann zeta function evaluated at argument n. The detailed phenomeno- 
logical consequences of the fully exponentiated 2- and 3-loop DGLAP-CS kernel set will 
appear elsewhere [19]. 
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3 Phenomenological Effects of IR-Improvement 



Let us now look into the effects of IR-improvement on the moments of the structure func- 
tions by discussing the corresponding effects on the moments of the parton distributions. 
In this section we work with the leading order results for definiteness of illustration. 

We know that moments of the kernels determine the exponents in the logarithmic 
variation [3-5] of the moments of the quark distributions and, thereby, of the moments of 
the structure functions themselves. To wit, in the non-singlet case, we have 

where ^ 

M^^t)= [ dzz'^-^'^^z^t) (50) 
Jo 



and the quantity A^^ is given (using (!39l) ) by 

Jo 



= CFFyps{%)e-^'''[B{n,^g) + B{n + 2,7,) - /,(7,)] (51) 
where B{x,y) is the beta function given by 

B{x,y) = T{x)T{y)/T{x + y) 
. This should be compared to the un-IR-improved result [3-5]: 



A"!'" = c, 



1 1 
-7: + 



n ■ 



2 n{n + 1) ^ j 



(52) 



The asymptotic behavior for large n is now very different, as the IR- improved expo- 
nent approaches a constant, a multiple of — /g, as we would expect as n — > 00 because 
lim^i^oo = for < 2 < 1 whereas, as it is well-known, the un-IR-improved result in 
(l52l) diverges as — 2Ci?lnn as n — > 00. The two results are also different at finite n: for 
n = 2 we get, for example, for as = .118 [23], 

j^NS ^ I Cf(-1.33) , un-IR-improved 
^ [C^(-0.966) , IR-improved 

so that the effects we have calculated are important for all n in general. For completeness, 
we note that the solution to 0491) is given by the standard methods as 



t,to large with t»to \ as{t) 
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where Ei{x) = dre'' /r is the exponential integral function, 

- ^^^Fypsiiy-^ [B{n, 7,) + B{n + 2, 7,) - fM] 



" /So 
a'. = a„, ( 1 



with 



^1 (ft.(^o) - «s(^)) 
2 ln(a,(to)/as(t)) 



(55) 



. We can compare with the un-IR- improved result in which the last line in eq. llMl) holds 
exactly with a'„, = 2A^^°/j3o. Phenomenologically, for n = 2, taking Qo = 2GeV and 
evolving to Q = lOOGeV, if we set Aqcd — ■2GeV and use Uf = 5 for definiteness of 
illustration, we see from eqs. fl54ll55p that we get a shift of the respective evolved NS 
moment by ~ 5%, which is of some interest in view of the expected HERA precision [24]. 

We stress that the size of the exponent 7^ is what one would expect from anal- 
ogy with QED [25], where with Q = 100 GeV we have the analogous result 7e = 
(^em/tt) (In Q^/mg — 1) = 0.054 whereas here, with = .118, so that it is about 10 
times Qem, we get a value for 7^ that is about 10 times gammae- There are no data 
that currently contradict the values we find for our 7^, = q, G. Evidently, from the exact 
result (IHj) we can consistently improve the values of the 7-,- as needed using more and 

more of the perturbative results for the functions SUMiji{QCD), Sqqd, Pn on the RHS 
accordingly, as we have noted. 

We give now the remaining elements of the leading anomalous dimension matrix in 
its 'best' IR-improved form for completeness: 



(56) 



A^^= [ dzz"-'P^''^iz)=2CGFypsilG)e'^'^{B{n + l,^G)+B{n + ^G-l,2) 
Jo 

+ ^{B{n + 1, + 2) + 5(n + 7G + 1, 2)) - /(7g)}, (57) 

2nfAf = 2nf [ dzz^~^P'J{z) = 2T{F)FYFs{lG)e^''' {B{n + 2, 1 + 7^) + B{n + 7^, 3)) 
Jo 

(58) 

where T{F) = ^Uf. We note that the un-exponentiated value of the last result in eq. fl58|) is 
a well-known one [3-5], 2T{F) ^(^^")"["^2) ' ^"^^ would be used whenever we do not choose 
to exponentiate Pqc- We will investigate the further implications of these IR-improved 
results for LHC physics elsewhere [19]. 
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4 Impact of IR-Improvement on the Standard Method- 
ology 



In sum, we have used exact re-arrangement of the QCD Feynman series to isolate and 
resum the leading IR contributions to the physical processes that generate the evolution 
kernels in DGLAP-CS theory. In this way, we have obviated the need to employ artificial 
mathematical regularization of the attendant IR singularities as the theory's higher or- 
der corrections naturally tame these singularities. The resulting IR-improved anomalous 
dimension matrix behaves more physically for large n and receives significant effects at 
finite n from the exponentiation. 

We in principle can make contact with the moment-space resummation results in 
Ref. [26] but we stress that these results have necessarily been obtained after making a 
Mellin transform of the mathematical artifact which we address in this paper. Thus, the 
results in Ref. [26] do not in any way contradict the analysis in this paper. 

We note that the program of improvement of the hadron cross section calculations for 
LHC physics advanced herein should be distinguished from the results in Refs. [10,27,28]. 
Indeed, recalling the standard hadron cross section formula 



where {Fi{x)} are the respective parton densities and a{xiX2s) is the respective reduced 
hard parton cross section, the resummation results in Refs. [27, 28] address, by summing 
the large logs in Mellin transform space, the X1X2 1 limit of (j{xiX2s) whereas the 
results above address the improvement, by resummation in x-space, of the calculation 
of the parton densities for all values of x. Thus, the program of improvement 

presented above is entirely complementary to that in Refs. [27, 28] and both programs 
of improvement are useful for precision LHC physics. The situation can be illustrated 
by comparing the results in Refs. [29] with our results herein. The key observation can 
already be made from eq.(2.1) in the latter paper, wherein it is made manifest that the 
resummation carried out therein, as an application of the methods in Refs. [27,28], is 
a resummation for the large N-Mellin space limit of the Mellin transform of the hard 
scattering coefficient function so that all of the IR effects in the parton densities are not 
included in this resummation. What we deal with here is however resummation of the 
IR effects in the kernels which generate exactly these IR effects in these parton densities 
directly in configuration space so that we work on a complementary aspect the formula 
(l59ll and this we do directly in x-space rather than in N-Mellin space. There is then 
no contradiction or repetition between our results and those in Ref. [29] . Refs. [27, 28] 
have stressed that the IR singularities of the parton densities are not addressed by their 
analysis. If we wanted to include the results in Refs. [27,28], as we have discussed in 
Refs. [15], we would make an inverse Mellin transform of their results and apply them to 
our resummation calculus for the hard cross section a (0:1X2 s) as it is illustrated in the 




(59) 
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Appendix; by construction, this would not affect the IR singularities for the evolution of 
the initial parton densities analyzed in this paper. 

Similarly, the analysis of the Pt distribution resummation in Ref. [10] deals with the 
hard scattering coefficient response to soft-gluon emission. The initial parton densities 
are input to this analysis. Evidently, as we work with the x-space initial parton density 
evolution in which any Pt has been integrated out, the analysis in Ref. [10] is entirely 
complementary to our work in this paper. 

Finally, we address the issue of the relationship between the re-arrangement that we 
have made of the exact leading-logs in the QCD perturbation theory and the usual treat- 
ment in the non-exponentiated DGLAP-CS theory. If one expands out the exponentiated 
kernels, using the distribution identity 



one can see that for example Pqq and P^^^ agree to leading order, so that the leading log 
series which they generate for the respective NS parton distributions also agree through 
leading order in —L where L is the respective big log in momentum-space. At higher 
orders then, we have a different result for the {Fi}, let us denote them by {F'i}, and a 
different result for the reduced cross section, let us denote it by a', such that we get the 
same perturbative QCD cross section. 



order by order in perturbation theory. The exponentiated kernels are used to factorize the 
mass singularities from the unfactorized reduced cross section and this generates a' instead 
of the usual a whose factorized form is generated using the usual DGLAP-CS kernels. For 
example, from ([7]), fl39|) and fl60l) . we see that, starting at order when one factorizes 

the mass singularities using P^^^ one makes a different subtraction from the respective 
fixed-order unfactorized hard cross section compared to what one would subtract if the 
factorization were done with the unexponentiated Pgg, resulting in a different reduced 
cross section, a'; but, the convloution of this different reduced cross section with the 
respective {-F/} that result from solving the DGLAP-CS equation with the IR-improved 
kernels would then give the same result for the hadron-hadron scattering cross section as 
one would get with the un- improved set {Fi, a}, order by order in perturbation theory. 

The entirely analogous statement holds for the structure functions in deep inelastic 
scattering. They are represented as a sum over products of the parton distributions and 
the attendant hard scattering cross sections. Again, if we use exponentiated kernels to 
generate the respective parton distributions, we use these kernels to isolate the respective 




(60) 




(61) 
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factorized hard scattering cross sections, giving the same perturbative QCD prediction 
order by order in perturbation theory for the respective deep inelastic structure functions. 

We thus have the same leading log series for a as does the usual calculation with 
un-exponentiated DGLAP-CS kernels. We have an important advantage: the lack of 
+-functions in the generation of the configuration space functions {-F'j, a'} means that 
these functions lend themselves more readily to Monte Carlo realization to arbitrarily soft 
radiative effects, both for the generation of the parton shower associated to the {F'i} and 
for the attendant remaining radiative effects in a' . 

We stress that, as on can see from our analysis above, unlike the standard MS un- 
exponentiated kernels, our exponentiated kernels contain powers of the product UgL, as 
these describe the large IR effects which we resum. These effects are then beyond the 
standard Wilson expansion, as they cause non-trivial modification of the moments of the 
exponentiated kernels relative to the moments of the MS kernels and it is the latter 
which correspond to the respective anomalous dimensions of the operators in Wilson's 
expansion when this expansion exists. We exchange the naive connection between the 
moments of the kernels and these anomalous dimensions for the improved IR behavior of 
the exponentiated kernels. We do this with an eye toward MC methods. 

The use of the new IR-improved kernels may seem as though one is over counting 
effects; after all, it is well-known [30] that the response of the DGLAP-CS equation for 
a delta-function at 2; = 1 for an initial value t = in the NS case is indeed a solution 
which exhibits exactly (1 — z)^'^'^ type behavior for the respective distribution that we 
have found for our NS exponentiated kernel. We are not over-counting this effect. For, 
consider the sequential application of the kernel P^^^{z) so that, first, a quark, splits 
via its action to a quark q{z) and gluons {G'j(^j)}, with Xlj'Ci = 1 — 2 in the standard 
notation, and, second, the latter quark, q{z), splits via the second application of P^^^{z') 
to a quark q{zz') and gluons {Gj(^j)}, with Ylij^'j = (1 ~ ^')^- the first action, the 
quantum transition was for g(l) — > q{z) on the quark Hilbert space while in the second 
it was for q{z) — »• q[zz'), a completely different quantum transition on the quark Hilbert 
space. The two transitions are completely independent in the leading log and using P^^^ 
to effect them does not over-count anything; it simply makes the description of each 
transition more accurate. 



5 Conclusions 

We have developed a new approach to precision QCD predictions for high energy colliding 
beam physics scenarios such as that afforded us by the impending turn-on of the LHC. 
We do not contradict any aspect of the traditional approach. We gain the advantages 
of improved IR behavior for the respective kernels, parton distributions and reduced 
cross sections, which should facilitate realization by multiple gluon MC methods. Such 
realizations for LHC physics will appear elsewhere. [19]. 
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Figure 2: The process Q'Q Q'" + Q" + n{G). The four-momenta are indicated in the 
standard manner: qi is the four-momentum of the incoming Q, q'2 is the four-momentum 
of the outgoing Q" , etc., and Q = u,d, s, c, b, G. 
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Appendix 

In this Appendix I, we present the new QCD exponentiation theory which has been 
developed in Refs. [15,16] as it is not generally familiar. The goal is to make the current 
paper self-contained. 

For definiteness, we will use the process in Fig. 2, Q' {pi)Q{qi) — > Q'" {p2)Q" {.(I2) + 
Gi{ki) ■ ■ ■ Gn{kn), as the proto-typical process, where we have written the kinematics 
as it is illustrated in the figure. This process, which dominates processes such as tt 
production at FNAL, contains all of the theoretical issues that we must face at the parton 
level to establish , as an extension of the original ideas of Yennie, Frautschi and Suura 
(YFS) [17], QCD soft exponentiation by MC methods - applicability to other related 
processes will be immediate. For reference, let us also note that, in what follows, we use 
the GPS conventions of Ref. [31] for spinors {u, v, u} and the attendant photon and gluon 
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polarization vectors that follow therefrom: 

with = and ( defined in Ref. [31], so that all phase information is strictly known 
in our amplitudes. This means that, although we shall use the older EEX realization 
of YFS MC exponentiation as defined in Ref. [32] , the realization of our results via the 
the newer CEEX realization of YFS exponentiation in Ref. [32] is also possible and is in 
progress [19]. 

Specifically, the authors in Refs. [16] have analyzed how in the special case of Born 
level color exchange one applies the YFS theory to QCD by extending the respective YFS 
IR singularity analysis to QCD to all orders in as- Here, unlike what was emphasized 
in Refs. [16], we focus on the YFS theory as a general re-arrangement of renormalized 
perturbation theory based on its IR behavior, just as the renormalization group is a general 
property of renormalized perturbation theory based on its UV(ultra-violet) behavior. We 
will thus keep our arguments entirely general from the outset, so that it will be immediate 
that our result applies to any renormalized perturbation theory in which the cross section 
under study is finite. 

Let the amplitude for the emission of n real gluons in our proto-typical subprocess, 
+ Q'°' Q"'^Q"'"' + n{G), where a, a, 7, and 7 are color indices, be represented by 

A<W"^ = ^m(5"^ (63) 

e 

is the contribution to TW*-"-* from Feynman diagrams with £ virtual loops. Sym- 
metrization yields 

^t'4/n p,),(,fi,.,,) Pr'^ (64) 

where this last equation defines symmetric function of its arguments arguments 

ki, ...,ki. A will be our infrared gluon regulator mass for IR singularities; n-dimensional 
regularization of the 't Hooft-Veltman [33] type is also possible as we shall see. 

We now define the virtual IR emission factor Sqcoik) for a gluon of 4-momentum k, 
for the A; — > regime of the respective 4-dimensional loop integration as in fl6^ . such that 



lira A; ip^^i°'{k)\ leading Casimir contribution SQCD{k)p%"'') = 0, (65) 

where we have now introduced the restriction to the leading color Casimir terms at one- 
loop0 so that in the expression for the respective one-loop correction p^'^'' and in that for 



^These correspond with maximally non-Abclian terms in Ref. [34] but computed exactly rather than 
in the eikonal approximation. 
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for Sqcoik) given in Refs. [16], only the terms proportional to Cp should be retained here 
as we focus on the ff— s-ff case, where f denotes a fermion. (Henceforth, when we refer to 
k gluons we are always referring for virtual gluons to the corresponding regime of 
the 4-dimensional loop integration in the computation of m]"\) 

In Ref. [16], the respective authors have calculated Sqcnik) using the running quark 
masses to regulate its collinear mass singularities, for example; n-dimensional regulariza- 
tion of the 't Hooft-Veltman type is also possible for these mass singularities and we will 
also illustrate this presently. 

We stress that SgcDik) has a freedom in it corresponding to the fact that any function 
ASqcoik) which has the property that \imk_,o k'^ AS QCD{k)pl)^'' = may be added to it. 

Since the virtual gluons in p^"'' are all on equal footing by the symmetry of this 
function, if we look at gluon i, for example, we may write , for k^ (0, 0, 0, 0) = O while 
the remaining ki are fixed away from O, the representation 

pP = SQcoike) * ■ ■ ■ , + /3l{ki, ■ ■ ■ , /c£_i; (66) 

where the residual amplitude /5^(fci, ■ ■ ■ , /c£_i; ki) will now be taken as defined by this last 
equation. It has two nice properties: 

• it is symmetric in its first i — 1 arguments 

• the IR singularities for gluon i that are contained in SqcDike) are no longer con- 
tained in it. 

We do not at this point discuss the extent to which there are any further remaining 
IR singularities for gluon i in Pl{ki, ■ ■ ■ , ki). In an Abelian gauge theory like QED, 
as has been shown by Yennie, Frautschi and Suura in Ref. [17], there would not be any 
further such singularities; for a non- Abelian gauge theory like QCD, this point requires 
further discussion and we will come back to this point presently. 

We rather now stress that if we apply the representation fl66p again we may write 

Pe''' = SQCD{h)SQCD{h^i) * p1"2(^1' " " " ' ^^-2) 
+SQCD{kt)f3l^i{ki, • • • , ki_2] ki-i) 
+SQCD{ki-i)(3}_i{ki, ■■■ , ki_2; h) 

+P]{ki, ■■■ , ki_2; ki-i, kg), (67) 

where this last equation serves to define the function /3|(/ci, • • • , /c£_2; fc^-i, kg). It has two nice 
properties: 

• it is symmetric in its first 1 — 2 arguments and in its last two arguments A;^_i, kg 

• the infrared singularities for gluons i — I and i that are contained in SQCD{ki~i) and 
SqcDike) are no longer contained in it. 



20 



Continuing in this way, with repeated apphcation of (I66p . we get finally the rigorous, exact 
rearrangement of the contributions to p^^^ as 

pf = SqcDih) ■ ■ ■ SQCD{h)(3l + E n SQCD{kj)P\{h) + ■■■ 

+f3f{h,--- ,k,), (68) 
where the virtual gluon residuals PKk'i, ■ ■ ■ ,k'^ have two nice properties: 

• they are symmetric functions of their arguments 

• they do not contain any of the IR singularities which are contained in the product 
SqcDik'i)--- Sqcd{K). 

Henceforth, we denote 13] as the function /5j for reasons of pedagogy. We can not stress 
too much that (j68p is an exact rearrangement of the contributions of the Feynman diagrams 

(n) 

which contribute to , it involves no approximations. Here also we note that the question 
of the absolute convergence of these Feynman diagrams from the standpoint of constructive 
field theory remains open as usual. Yennie, Frautschi and Suura [17] have already stressed that 
Feynman diagrammatic perturbation theory is non-rigorous from this standpoint. What we do 
claim is that the relationship between the YFS expansion and the usual perturbative Feynman 
diagrammatic expansion is itself rigorous even though neither of the two expansions themselves 
is rigorous. 

Introducing ()68p into (I63p yields a representation similar to that of YFS, and we will call it 
a "YFS representation", 

oo 

J^[n) ^ ^asBQCD (69) 

i=o 

where we have defined 

/d'^k 
■(^Zrj2-^^QCDik) (70) 

and 

We say that (f69]) is similar to the respective result of Yennie, Frautschi and Suura in Ref. [17] and 
is not identical to it because we have not proved that the functions f3i{ki, ki) are completely 
free of virtual IR singularities. What have shown is that they do not contain the IR singularities 
in the product SqcDiki) ■ • • SQCDih) so that m^"'^ does not contain the virtual IR divergences 
generated by this product when it is integrated over the respective 4j-dimensional j-virtual gluon 
phase space. In an Abelian gauge theory, there are no other possible virtual IR divergences; in 
the non-Abelian gauge theory that we treat here, such additional IR divergences are possible 
and are expected; but, the result ([69]) does have an improved IR divergence structure over ([63]) 
in that all of the IR singularities associated with Sqcoik) are explicitly removed from the sum 
over the virtual IR improved loop contributions m^"^ to all orders in as{Q)- 
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Turning now to the analogous rearrangement of the real IR singularities in the differential 
cross section associated with the A4^^\ we first note that we may write this cross section as 
follows according to the standard methods 

= n\ ./ ii ThTW^ 



m=l ^ ™ ' " ' i=l 

'{Pi,qi,P2,q2,ki, - ■ ■ ,kn) fTn"' ('2) 



where we have defined 



oo 

P^''Hpi,qi,P2,q2,hr-- ,kn)= (^3) 

color, spin j=0 

in the incoming QQ' cms system and we have absorbed the remaining kinematical factors 
for the initial state flux, spin and color averages into the normalization of the amplitudes 
A<(") for reasons of pedagogy so that the p*^"^ are averaged over initial spins and colors 
and summed over final spins and colors. We now proceed in complete analogy with the 
discussion of p^"^ above. 

Specifically, for the functions p^'^\pi, qi,P2, Q2, ki, - ■ ■ , kn) = p^"'\ki, ■ ■ ■ , kn) which are 
symmetric functions of their arguments /ci, ■ ■ ■ , kn, we define first, for n = 1, 

lim P ( p'-^'^H^) I leading Casimir contribution — 'S'qc'D (^)P*-°M = 0, (74) 

where the real infrared function Sqcnik) is rigorously defined by this last equation and is 
explicitly computed in Refs. [16], wherein we retain here only the terms proportional to 
Cp from the result in Ref. [16] ; like its virtual counterpart SqcDik) it has a freedom in 
it in that any function ^Sqcnik) with the property that lim|g|^Q A;^AS'qc'£)(/c) = may 
be added to it without affecting the defining relation (174|) . 

We can again repeat the analogous arguments of Ref. [17], following the corresponding 
steps in fl66|) - fl7T]) above for Sqcd to get the "YFS-like" result 



oo 



_ SUMir{QCD)V- lf\'^h. f J^-(P^+cii-P9-ci9-y.k,)+Doc.T^ 

■ Pn{kl, . . . , kji) rj rj- 

P2 92 



with 



SUMiBiQCD) = 2a,ReBQCD + 2asBQCD{Km..) 

/(fik ~ 
-j^Sq^CT>{k)e{K^s^ - k), 
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—SqcBik) \e-'y^ - e{K^^^ - k)\ , (76) 
^/3o = da(i-i°°P) - 2asReBQCDdaB, 

^Pi = da''^ -SQCDik)daB, ... (77) 

where the /3„ are the QCD hard gluon residuals defined above; they are the non-Abehan analogs 
of the hard photon residuals defined by YFS. Here, for illustration, we have recorded the relation- 
ship between the n = 0, 1 through 0{as) and the exact one-loop and single bremsstrahlung 
cross sections, da^^~^°"'^\ da^^ , respectively, where the latter may be taken from Ref. [35] We 
stress two things about the right-hand side of ([75]) : 

• It does not depend on the dummy parameter K^ax which has been introduced for can- 
cellation of the infrared divergences in SU Mir{QCD) to all orders in as{Q) where Q is 
the hard scale in the parton scattering process under study here. 

• Its analog can also be derived in our new CEEX [32] format. 

We now return to the property of (|75p that distinguishes it from the Abelian result derived 
by Yennie, Frautschi and Suura - namely, the fact that, owing to its non-Abelian gauge theory 
origins, it is in general expected that there are infrared divergences in the /3„ which were not 
removed into the Sqcd, Sqcd when these infrared functions were isolated in our derivation of 

m. 



More precisely, the left-hand side of (1751) is the fundamental reduced parton cross 
section and it should be infrared finite or else the entire QCD parton model has to be 
abandoned. 

There is an observation in the literature [36] that unless we use the approximation 
of massless incoming quarks, the reduced parton cross section on the left-hand side of 
flTSl) diverges in the infrared regime at 0{al{Q)). We do not go into this issue here but 
either use the quark masses strictly as coUinear limit regulators so that they are set to 
zero in the numerators of all Feynman diagrams in such a way that the limit lim^2/£;2^o; 
where Eg is the quark energy, is taken everywhere that it is finite or, alternatively, we 
use n-dimensional methods to regulate such divergences while setting the quark masses to 
zero as that is an excellent approximation for the light quarks at FNAL and LHC energies 
- we take this issue up elsewhere. 

From the infrared finiteness of the left-hand side of fl75p and the infrared finiteness of 
SUMib{QCD), it follows that the quantity 

d^ = p-SUMiR(QCD) 

must also be infrared finite to all orders in a^. 

As we assume the QCD theory makes sense in some neighborhood of the origin for 
as, we conclude that each order in must make an infrared finite contribution to ciiJexp. 
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At 0{a'^{Q)) , the only contribution to (io"exp is the respective Born cross section given 
by /^o'^'' in fl75l) and it is obviously infrared finite, where we use henceforth the notation 
to denote the part of Thus, we conclude that the lowest hard gluon 

residual f3^^ is infrared finite. 

Let us now define the left-over non-Abelian infrared divergence part of each contribu- 
tion l3n^ via 

where the new function jsli' is now completely free of any infrared divergences and the 
function D/3n^ contains all left-over infrared divergences in (3n^ which are of non-Abelian 
origin and is normalized to vanish in the Abelian limit fate where fabc are the group 
structure constants. 

Further, we define D/3n^ by a minimal subtraction of the respective IR divergences 
in it so that it only contains the actual pole and transcendental constants, 1/e — Ce for 
6 = 2 — d/2, where d is the dimension of space-time, in dimensional regularization or In 
in the gluon mass regularization. Here, Ce is Euler's constant. 

For definiteness, we write this out explicitly as follows: 



dPh DPi'^ = Y,dfln\\' 



where the coefficient functions d"'^ are independent of A for A ^ and dPh is the 
respective n-gluon Lorentz invariant phase space. 

At 0{a^{Q)), the IR finiteness of the contribution to ddcxp then requires the contri- 
bution 

- / E ? n i ^s,oo(k,) E i n 

e=0 ' j = l-^l'j>^ni,ax 3 j = ' j=e+i 
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to be finite. 

From this it follows that 



e=Q ■ j=l -'^J^^maa; J j^Q ' j=e+l 

d^^J r,a(n-e-i).,^ ,^ .d^P2d^q2 



>ike+l^■■■^h+^)^^ (79) 



P2 Q2 



is finite. Since the integration region for the final particles is arbitrary, the independent 
powers of the IR regulator In(A^) in this last equation must give vanishing contributions. 
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This means that we can drop the D(3n from our result (1751) because they do not make a 
net contribution to the final parton cross section ^exp- We thus finally arrive at the new 
rigorous result 



d^exp = ^ da"' 



„SUMir(QCD) 



E/n 



A k, J (27r) 



(80) 



where now the hard gluon residuals f3n{ki , ■ ■ ■ , kn) defined by 



f^nikl ^ ■ ■ ■ ikn — ^ ^ (in (^1 ) • • • ) kn) 
t=0 



(81) 



are free of all infrared divergences to all orders in as{Q)- This is a basic result of this Appendix. 

We note here that, contrary to what was claimed in the Appendix of the first paper in 
Refs. [16] and consistent with what is explained in the third reference in [16], the arguments in 
the first paper in Refs. [16] are not sufficient to derive the respective analog of eq. (l80p : for, they 
did not really expose the compensation between the left over genuine non-Abelian IR virtual and 
real singularities between / dPhPn and f dPhPn+i respectively that really distinguishes QCD 
from QED, where no such compensation occurs in the (5n residuals for QED. 

We point-out that the general non-Abelian exponentiation of the eikonal cross sections in 
QCD has been proven formally in Ref. [34]. The contact between Ref. [34] and our result (j80p 
is that, in the language of Ref. [34], our exponential factor corresponds to the N=l term in the 
exponent of eq.(lO) of the latter reference. One also sees immediately the fundamental difference 
between what we derive in (I80p and the eikonal formula in Ref. [34]: our result (I80p is an exact 
re-arrangement of the complete cross section whereas the result in eq.(lO) of Ref. [34] is an 
approximation to the complete cross section in which all terms that could not be eikonalized 
and exponentiated have been dropped. 
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